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Abstract: We consider non-universal 'minimal' Z' models, whose additional U(l) 
charge is a non- anomalous linear combination of the weak hypercharge Y , the baryon 
number B and the partial lepton numbers (L e , L M , L r ), with no exotic fermions be- 
yond three standard families with right-handed neutrinos. We show that the observed 
pattern of neutrino masses and mixing can be fully reproduced by a gauge-invariant 
renormalizable Lagrangian, and flavor-changing neutral currents in the charged lep- 
ton sector are suppressed by a GIM mechanism. We then discuss the phenomenology 
of some benchmark models. The electrophilic B — 3L e model is significantly con- 
strained by electroweak precision tests, but still allows to fit the hint of an excess 
observed by CDF in dielectrons but not in dimuons. The muonphilic B — 3L M model 
is very mildly constrained by electroweak precision tests, so that even the very early 
phase of the LHC can explore significant areas of parameter space. We also discuss 
the hadrophobic L M — L T model, which has recently attracted interest in connection 
with some puzzling features of cosmic ray spectra. 
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1. Introduction 

Extensions of the Standard Model (SM) at the TeV scale, with an additional U(l) 
factor in the gauge group associated with a heavy neutral gauge boson Z' ', have 
often been considered, with various theoretical motivations (for recent reviews and 
references, see e.g. [1]). 

In a recent paper [2], three of us discussed the phenomenology of minimal Z f 
models, previously introduced in [3] and identified as the most economical U(l) 
extensions of the SM that do not spoil renormalizability. Making reference to the 
SM particle content, the key ingredients of minimality are: no exotic vectors, apart 
from a single Z f associated with a U(l) factor in the gauge group, commuting with 
Gsm = SU(3)c x SU(2) L x U(1)y] no exotic fermions, apart from one right-handed 
neutrino, singlet under Gsm, for each of the three SM families; anomaly cancellation. 

An additional simplifying assumption made in [2] was the family-independence 
of the additional U(l) charge, which could then be identified with an arbitrary linear 
combination of the weak hypercharge Y and B — L, where L = L e + L M + L r is the 
total lepton number. In such a case, electroweak precision tests (EWPT) strongly 
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constrain the parameter space of the models, leaving little virgin land to be explored 
in the very early phase of the LHC, when energy and luminosity will be limited [4] 
with respect to the design parameters. 

In this paper, we relax the assumption made in [2], allowing the additional 
U(l) generator to be an arbitrary non-anomalous linear combination of the weak 
hypercharge Y ) the baryon number B and the partial lepton numbers L e , L M , L T . It 
is immediate to check (see, e.g., [5]) that, when baryon and lepton numbers appear 
only in a linear combination X of the three generators (B — 3L a ) (a = e,/x, r), 
any such model is anomaly-free. Indeed, it is easy to check that, if we add the 
requirement that all fermion masses and mixing be generated by a renormalizable 
gauge-invariant Lagrangian, with a suitable Higgs content, then this is the most 
general solution. Models have already been considered with gauged B — 3L T [6], 
B — 3L M [7] and L a — (a ^ b) [9]. In particular, models with gauged B — 3L M were 
suggested 1 as possible non-SM explanations [7] of the NuTeV 'anomaly' [10], and 
models with gauged L a — (a ^ b) as non-SM explanations [11] of possible excesses 
in e ± cosmic ray spectra [12]. Some constraints on the above models from EWPT 
were discussed in [13], and other anomaly free models, where neutrino masses are 
generated by non-renormalizable gauge-invariant interactions, were considered in [5]. 

Our paper is organized as follows. In Section 2, we discuss the theoretical in- 
put for the following phenomenological analysis. We first review the structure of 
the fermionic neutral currents. We then show that the observed pattern of fermion 
masses and mixings, in particular the one in the neutrino sector, can be obtained 
from a gauge-invariant renormalizable Lagrangian, with no more fine-tuning of the 
Yukawa couplings than the one already needed in the SM to reproduce the observed 
electron mass. Moreover, all flavor-changing neutral current processes (FCNC) in- 
volving charged leptons are strongly suppressed by the neutrino mass differences, 
as in the Glashow-Iliopoulos-Maiani (GIM) mechanism [14] of the SM. We end this 
theoretical section by studying the renormalization group equations (RGE) for the 
effective couplings controlling the weak neutral currents, as well as the Z-Z f mix- 
ing. Assuming that the model is valid up to very high scales, such as some super- 
unification or grand-unification scale, we identify a favored region for the effective 
low-energy couplings and mixing. In Section 3 we review the general aspects of the 
phenomenology of the class of models under consideration. We start with the con- 
straints from electroweak precision tests, where we complete and update the fit of 
[15], based in turn on [16] and used in [2] for constraining the universal models where 
X = B — L. We then briefly recall the procedure followed in [2] for extracting the 
Tevatron bounds and assessing the discovery prospects for the very early phase of 
the LHC, to be repeated here for non-universal minimal Z 1 models. In Section 4 we 
apply our results to the study of three benchmark models that we find particularly 

1 One of us (A.S.) thanks C. Boehm for discussions about why, after correcting a sign mistake, 
B - 3L M does not fit the NuTeV anomaly [8]. 
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interesting: 

1. an electrophilic model, corresponding to X = B — 3L e , which could explain, 
for values of the parameters not excluded by EWPT, a small excess observed 
at the Tevatron in the dielectron sample [17, 18], at invariant masses around 
240 GeV, but not in the dimuon sample [19]; 

2. a muonphilic model, corresponding to X = B — 3L M , which is subject to much 
milder constraints from EWPT, thus has a much wider area of parameter space 
accessible, via the dimuon signature, already in the very early phase of the LHC. 
In a sense, this model is another example of the supermodels recently discussed 
in [20], defined as those for which the LHC sensitivity with only 10 pb _1 is 
greater than that of the Tevatron with 10 fb -1 ; 

3. a hadrophobic model, corresponding to X = L M — L r , which has recently 
attracted considerable interest [11] as a possible non-SM explanation of the 
positron excess in cosmic ray data [12]. 

Finally, in Sect. 5 we present our conclusions. 



2. Theory 

2.1 Parameterization 

Along the lines of [2], we work in a field basis where gauge boson kinetic terms are 
canonical and gauge boson masses are diagonal, and write the two fermionic currents 
coupled to the neutral massive gauge bosons as: 

r z = cos & J£ - sin 0' J£ /0 , J£, = sin 0' + cos & j£ /0 , (2.1) 

where 

J" z o =9zJ2lr (T 3L - sin 2 d w Q) f, (g z = , (2.2) 

/ 

is the SM expression for the current coupled to the SM Z° (we recall that, in the 
presence of mixing, Z° does not coincide with the mass eigenstate Z), and 

Jz>o = ^^(9yY + 9x X) f = gzY,f'fQ*f- ( 2 - 3 ) 
/ / 

In the above expressions, / = {u La) d La) u Ra) d Ra) v La) e La) v Ra) e Ra ) runs over the 
different chiral projections of three families of SM fermions, Q = T 3L + Y is the 
electric charge and 

X= ]T ^-(B-3L a ) = pB-X e L e -X^-X T L T , /?=I^A a (2.4) 
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Table 1: The charges of the SM fermions controlling the weak neutral currents. 



where the A a are three arbitrary real coefficients. The charges of the SM fermions 
needed for evaluating the currents of Eqs. (2.1)-(2.3) are collected 2 in Tab. 1, where 
for expressing Q z > we make use of the ratios 

~ _ 9y ~ _ 9x , n rA 

9y = — , 9x = — • (2.5) 
9z 9z 



Since the SM Higgs doublet H has 3 X = 0, and for breaking X we consider 
only S r £/(2) j £ / -singlet Higgs fields, which must have Y = Q — T 3L = 0, to avoid the 
spontaneous breaking of the electric charge, we can express the Z-Z f mixing angle 9 f 
in terms of gy and M z >, 

M 2 n 2 n 2 

tmd> = -gy f f where Mf = ^ (2.6) 



2 Notice that, with respect to [2], we switch from a basis of left-handed fermions and antifermions 
to a basis of left- and right-handed fermions. 

3 This property is shared by the MSSM Higgs doublets Hi and H2. 
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is the SM expression for the Z° mass. Notice that 6' ^ whenever g Y ^ 0, because 
in such a case gauge invariance of the SM Yukawa terms forces the SM Higgs to be 
charged under the extra £7(1), thus producing a Z-Z' mixing. 

We can then study the Z' phenomenology in terms of five unknown parameters: 

• the Z' mass M z > and the two coupling constants (gy^gx) or, equivalently, 

• two independent ratios of the three flavor A a coefficients (the overall scale can 
be reabsorbed in g x ) that define X in Eq. (2.4). 

We will mention later possible additional parameters describing the enlarged Higgs 
sector and the right-handed neutrino masses, but we will mostly consider physical 
contexts in which these additional parameters play a negligible role. 

2.2 Fermion masses and mixing 

To make the model realistic, we have to make sure that the observed pattern of 
fermion masses and mixing can be generated. Our philosophy is to stick to renormal- 
izable interactions, formally sending the mass scale of non-renormalizable operators 
to infinity, and to introduce the minimal number of Higgs fields that can realize the 
required gauge and flavor breaking. 

We gauge baryon number B rather than baryon flavors, so that, in the quark 
sector, masses and mixing can be generated exactly as in the SM, by 

- £g k = u^Y u q L H + T R Y D q L H + h.c. , (2.7) 
where obvious contractions of SU(3), SU{2) and family indices are understood, 

H=( H *o) -(+1A0), H = to*rf = (_^°_) ~ (-1/2,0) (2.8) 

are the SM Higgs field and its conjugate (the numbers in brackets are the charges Y 
and X ) respectively), Y u and Y D are 3x3 complex matrices, unrestricted by the 
additional U(l) gauge invariance. As in the SM, the transition from the interaction 
basis to the mass basis for the quark fields generates the CKM matrix but preserves 
the diagonal form of the quark contribution to all tree-level gauge neutral currents. 

We then write down the most general Yukawa couplings involving the SM Higgs 
in the leptonic sector. Their general form will be as in the SM, 

- £g, k = e^Y E l L H + 1^Y N l L H + h.c. , (2.9) 

but, since we gauge lepton flavor rather than lepton number, some of the entries in 
the complex 3x3 matrices Y E and Y N can be now forbidden by the extra £7(1). To 
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understand this point, we can treat Y E and Y N as spurions, and write down the X 
charges of their individual entries: 

X(Y a ^=X(Y a ^) = \ b -X a . (2.10) 

Therefore, diagonal entries are always allowed, whilst those off-diagonal entries for 
which A5 7^ X a are forbidden. Anyway, there is no obstruction to the generation of 
the observed values of the charged lepton masses. 

For neutrino masses, we choose to make use of the (type-I) see-saw mechanism, 
to avoid fine-tuning the Yukawa couplings in Y N to a higher extent than those in 
Y E . We then allow for right-handed neutrino mass terms of the form 

Cm = \¥r)M r (v) V^ t + h.c. , (2.11) 

where M R (ip) is a complex symmetric 3x3 matrix that may contain bare mass terms 
as well as terms proportional to non-SM Higgs fields, singlets under the SM gauge 
group but carrying a non-vanishing X charge. As before, we can treat the entries of 
Mr(<p) as spurions and write: 

X[M R (i P ) ab ] = X a + X b . (2.12) 

For entries with X = 0, bare (field-independent) neutrino mass terms are allowed, 
whose size is not restricted by gauge invariance and unitarity. For entries with 1^0, 
mass terms can be generated by introducing renormalizable Yukawa couplings to the 
Higgs (es) (fx that break the extra U(l), with charge X under the extra U(l) and 
neutral under the SM gauge group: 

¥*~(0,X), (2.13) 

in the same notation of Eq. (2.8). Therefore, depending on their X charges, some 
right-handed neutrinos might acquire mass only after the breaking of the extra U(l). 
In conclusion, the right-handed neutrino masses can be of order M z > or higher. The 
three heavy neutrinos can be made sufficiently heavy to become negligible in the 
discussion of the Z' phenomenology. Similarly, the physical degrees of freedom in 
the non-SM Higgs fields of Eq. (2.13), whose vacuum expectation values (vevs) give 
mass to the Z' and to right-handed neutrinos, can be made sufficiently heavy and/or 
sufficiently decoupled from the SM Higgs field to be also negligible in the present 
discussion of the Z' phenomenology. 

We conclude this section by showing that our minimal non-universal renormal- 
izable 4 Z' models are compatible with the observed neutrino masses and mixing 

4 Here we differ from [5], which aimed at forbidding neutrino mass terms from renormalizable 
couplings, in the attempt of explaining in such a way the smallness of neutrino masses. We insist 
instead, coherently with the requirement of anomaly cancellation, on the renormalizabiliy of the 
theory at the TeV scale, meaning that all non-renormalizable operators allowed by the symmetries 
are suppressed by a very high physical cut-off scale. We do not see a problem, in the present context, 
in having the Yukawa couplings in Y N of the same of order of the electron Yukawa coupling in Y E . 
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(for a review, see, e.g., [21]), which turn out to contain, as in the SM, all lepton- 
flavor violation of experimental interest: the new Z 1 gauge interactions do not lead 
to dangerous FCNC processes involving charged leptons, which get suppressed by 
a GIM-like mechanism [14]. We want to reach the usual basis where the charged 
lepton mass matrix is diagonal and the mass matrix for the three light neutrinos is 
written in terms of the mass eigenvalues mi ;2 ,3 and the mixing matrix £7, as 

m v = U* • diag(mi, ra 2 , m 3 ) • . (2-14) 

So, we first diagonalize the charged lepton mass matrix, 

M E — Y E (H°) , (2.15) 

by performing the usual unitary transformations 

l L ^U L l L) e R ^U R e R) (2.16) 

where both Ul and U R are 3x3 unitary matrices, so that 

U ] R M E U L = diag(m e , m M , m T ) , (2.17) 

and we can identify the redefined v L as the neutrino flavor eigenstates, partners of the 
charged lepton mass and flavor eigenstates in SU(2) L doublets. Correspondingly, the 
Yukawa couplings associated to Dirac neutrino masses get redefined as Y N ^Y N Ul 
After the rotation the coupling of the Z' to the leptons reads 

9zZ'^ (T Ll »UiQ z ,U L l L + ^U R Q z ,U R e R + V^fQz^ R ) • (2.18) 

The important observation is that Ul and Ur do not mix sectors with different X 
charges, so that flavor mixing is present only between flavors with the same X charge, 
thus with the same Z' interaction, and no tree-level FCNC involving the SM charged 
leptons are generated. Flavour changing neutral currents do appear in the neutrino 
sector, as a consequence of the rotation required for diagonalizing the neutrino mass 
matrix. Such FCNC are however suppressed by the light neutrino masses, as in the 
SM, and do not produce dangerous effects. 

We now show that this restricted framework, leaving generic the right-handed 
neutrino mass matrix Mr, has enough freedom to reproduce the observed light neu- 
trino masses and mixing. The standard see-saw formula 

rrf = (M N ) T • M R X • M N , (2.19) 

where M N = Y N (H°), can be trivially solved for M R = (M N ) T • (m^)" 1 • (M N ), 
when det m v ^ 0, meaning that any 5 neutrino mass matrix m v can be obtained from 
any given Dirac mass matrix M N for an appropriate M R . 

5 The argument remains of course true in the case of one massless neutrino, still allowed by the 
data. 
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2.3 Renormalization group effects 

We recall that we denoted as g and g' the SM gauge couplings of SU(2) L x U{1) Y) 
and as gx and gy the gauge couplings of the X and Y components of the extra 
U(l) current, see Eq. (2.3). A given U(l) is specified by assigning X, see Eq. (2.4). 
The values of these couplings are scale-dependent, and their scale-dependence is 
controlled by the corresponding RGE. In the formalism of [2], but referring this time 
directly to the three U(l) gauge coupling constants, we can write at one-loop: 

dg x 1 



dt " 16tt 2 



(b X x gx + 2 b Y x 9x 9y + b YY gx 9y) , (2.20) 

% = —^ YY (j\ (2.21) 
dt 16tt 2 YYy ' y J 

^ = vL^ ( brY 9y ^ + 29 ^ + 2 brx 9x ^ + 9 ^ + bxx 9x 9y ) ' (2,22) 

where t = log(Q/Qo)? Qo is a reference scale, 

= Of Of + \Y1 Of Of . (A 5 = y, x) , (2.23) 

/ s 

and / and s are the two-component fermions and the complex scalars in the theory, 
respectively. 

In the following, when discussing our benchmark models, we will identify some 
GUT-favored regions in the (gy^9x) plane according to the following procedure. 
For choosing the boundary conditions at the GUT scale M U) we normalize all U(l) 
charges as T 3L , over three fermion generations with right-handed neutrinos, and we 
take Mjj = 10 16 GeV as a reference value. In typical GUTs, Mjj can vary within 
approximately two decades around such reference value, but the difference in our 
estimate of the GUT-favored region is of the order of other threshold effects that 
we reabsorb in the wide ranges we assume below for other parameters. Then, we 
compute the boundary value g'(Mu) using the phenomenological input g'(M z ) = 
e(M z )/cos6 w (M z ), with a~^(M z ) ~ 128 and sin 2 6 W (M Z ) ~ 0.23, and the SM one- 
loop RGE. We then allow the Z f coupling at the unification scale au = g^/(47r) = 
g\,(Mij) I (Atk\ to vary within the generous range 

Taking into account that the SM RGE would predict au ~ 1/45, our upper and 
lower bounds leave a margin of more than a factor of two to account for threshold 
corrections, new particles at the TeV scale and other model-dependent effects. Cor- 
respondingly, we determine the GUT-favored region of the (gy^dx) plane by making 
use of the one-loop RGE of eqs. (2.20)-(2.23): when discussing different benchmark 
models, the result will be presented as a colored band. In addition, we will study the 
stability with respect to the RGE evolution of the C pure-X' models, in which g Y = 
at the unification scale M v . 
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3. Phenomenology 



3.1 Electroweak precision tests 

As already mentioned, EWPT can be used to constrain the parameter space of 
minimal Z* models. The analysis of [2] was limited to the universal case where 
X = B — L, and was done under some approximations that allowed to make use of 
the results of [15] without performing a new fit. In particular: instead of the full set 
of relevant electroweak precision data, only nine pseudo-observables were used; the 
experimental and theoretical inputs of the fit were not updated beyond Winter 2006; 
the constraints coming from the muon anomalous magnetic moment (for a review, 
see, e.g., [22]) and from the NuTeV experiment [10] on neutrino-hadron deep inelastic 
scattering were not included. 

Here we perform a new, full global fit, employing the most recent measure- 
ments of the following precision observables performed at LEP1, Tevatron, SLC 
and other facilities: pole masses M Z) M w = (80.399 ± 0.023) GeV [23] and m t = 
(173.1 ± 1.3) GeV [24]; Fermi constant for fi decay; strong and electromagnetic cou- 
pling at M Z ] total Z width; e + e~ hadronic cross section at the Z peak, forward- 
backward asymmetry in the £ (= e, /x, r), 6, c final states; r polarization asymmetry; 
Z branching fraction into hadrons, bb ) cc\ left-right polarization asymmetry in the £ 
(= e, /x, r), 6, c final states. 

We also include the following LEP2 data, measured at various values of ^fs be- 
tween 183 and 207 GeV: a(e + e~ — > qq)] cross-sections and forward-backward asym- 
metries for e + e~ — ► t + t~, bb] angular distribution do(e r e~ —> e + e~)/dcos# 
as measured by OPAL and ALEPH. 

We additionally include the following low-energy measurements, which are less 
relevant in the fit to universal models, but can play a role in some region of the 
parameter space of non-universal models: M0ller scattering at Q 2 = 0.026 GeV 2 ; 
atomic parity violation (APV) in Cs as re-analyzed in [25]: Qw = —73.16 ±0.35; gL 
and qr in neutrino-nucleon scattering as measured by NuTeV; anomalous magnetic 
moment of the muon [22]. The last three observables are plagued by controversial 
theoretical uncertainties: we adopt the standard estimates and will warn the reader 
whenever these less relevant measurements play a significant role in the global fit. 

In the above lists, we provided numerical values and references only for those 
observables that have been updated with respect to [15, 16]. 

The corrections to all the above observables are computed by integrating out the 
heavy gauge boson coupled to the total (fermionic and Higgs) J^ /0 current, obtaining 
the effective Lagrangian: 

£ eff -£sM ' (3 ' 1} 
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Regions allowed by EWPT Regions allowed by EWPT 




Figure 1: Universal model, X=B—L. Left: The regions allowed by all EWPT at 95% CL, 
for M z , = 200, 500, 1000, 1500, 2000, 2500 and 3000 GeV (from inner to outer). Right: The 
regions allowed at 95% CL, for M z > = 1 TeV, by LEP2 (yellow), LEPl-APV-NuTeV-SLD 
(blue) and all combined (red). 

valid under the assumption that we are perturbing around the SM fit. The term pro- 
portional to \H^ iD )xH\ 2 1 'M\, affects the Z mass, whilst the squared fermion current 
gives four-fermions operators. The mixed term affects the Z couplings to fermions. 

The results of our improved and updated fit, applied to the parameter space of 
the minimal universal Z 1 model with X = B — L, discussed in [2], are displayed in 
Fig. 1: for some representative values of the Z f mass, the corresponding 95% CL 
allowed regions in the (gy^9x) plane are shown. On the left-hand side, the region 
allowed by the global fit is the red one enclosed by each contour, and the different 
mass values are (from inner to outer): M z , = 200, 500, 1000, 1500, 2000, 2500 and 
3000 GeV. Dashed lines and empty (full) dots remind us of the GUT-favored region 
and of the GUT-inspired (-derived) benchmark models, whose precise definitions can 
be found in [2]. 

Comparing with the corresponding figure of [2], we see that our improved and 
updated fit gives slightly more stringent bounds, and that the ellipsoidal curves de- 
limiting the allowed regions are now less tilted with respect to the vertical axis. A 
detailed check shows that the effect is mostly due to the transition from an approxi- 
mate to a full fit, and is in agreement with the findings of [15]. The updates of some 
input data such as m t and M W) and the inclusion of some new ones such (g — 2) M 
and NuTeV, do not affect the result in a significant way. To understand better the 
different effects, the right-hand side of the figure shows, for the representative value 
M z > = 1 TeV, how the regions allowed by LEP2 (yellow) and LEPl-APV-NuTeV- 
SLD (LEP1+ for short, blue) combine to give the globally allowed (red) region. 
The behavior is the expected one: data at the Z peak and at lower energies 6 are 
mostly sensitive to mixing effects, thus the blue band surrounds the vertical axis 

6 We are dealing here with a vector-like charge X for the extra U(l). 
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Table 2: 95% CL bounds (1-parameter fit) on the Z f masses from EWPT, corresponding 
to the specific universal models {X — B — L) represented by the nine points in Fig. 1. The 
^/(°) mo d e i s are those represented by empty points, while Z'^ (Z f<yW ^) corresponds to the 
three external (internal) black points; see [2] for details. 



that corresponds to vanishing mixing; data at LEP2 energies are mostly sensitive to 
Z'-exchange diagrams with some external electrons, thus the axis of the yellow band 
roughly corresponds to the direction that minimizes the Z' couplings to electrons. 

As in [2], we report in Tab. 2 the bounds on the Z 1 masses (in a 1-parameter 
fit) for some particular choices of the couplings corresponding to GUT-inspired and 
SUSY-GUT benchmark models. We observe that, because of the slightly different 
inclination of the ellipsoidal curves with respect to those from the approximate fit 
of [2], the bounds on the Q B — L models' are now slightly less stringent, whilst the 
bounds on the c x _m °dels' and on the c T 3jR -models' are now slightly more stringent. 

The validity of our fit to EWPT extends to the non-universal minimal models 
considered in this paper: the corresponding bounds on the parameter space will be 
illustrated in Sect. 4 for the three chosen benchmark models. 

Before concluding this discussion, we would like to point out that the bounds 
obtained in our fit are consistent with those given in [26], but considerably more 
stringent than those obtained by other authors (see, e.g., [1]). We identify two 
specific reasons for this difference. First, some authors do not take the constraints 
from LEP2 into full account; the right-hand side of Fig. 1 shows that LEP2 data are 
important. Second, some authors treat the Z-Z' mixing angle 9 f as an independent 
free parameter, the ^-minimization then selects a special value of 6 f very close to 0, 
which strongly relaxes the bounds from EWPT. However, for any given model, the 
value of 9 f is fixed by the gauge couplings and by the Higgs vevs to a value that, 
in general, is different from zero: in our class of models, it is fixed by Eq. (2.6). In 
models with a more exotic Higgs content, the cancellation of the mixing angle can 
be achieved only by fine-tuning the gauge couplings and the Higgs vevs and would 
be scale-dependent. Without such fine-tuning, also in those models the Z' bounds 
from EWPT would become much stronger, of the order of those reported in Tab. 2. 

3.2 Direct searches at hadron colliders 

We briefly recall here the general procedure followed in [2] for extracting the bounds 
from direct searches at the Tevatron and for assessing the discovery prospects in the 
very early phase of the LHC, to be repeated here for non-universal minimal Z' models. 
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This will prepare the ground for Section 4, where we will apply such a procedure to 
the study of our three benchmark models and comment on their peculiar features. 

Both at the Tevatron and at the LHC, we concentrate on the processes pp (pp) — > 
(Z 1 —> £ + £~) + X, (I = e, fi): these are two very clean channels to look for, and after 
simple generous cuts the irreducible background is dominated by the well-understood 
SM Drell-Yan (DY) processes. For the class of models under consideration, other 
channels into SM final states are not expected to be competitive for exclusion or 
discovery, even if they could play a role in the determination of the Z' couplings 
after a future discovery. 

For any given model in our class, specified by the choice of X, we compute 
the Tevatron Z 1 production cross-section multiplied by the branching ratio into two 
charged leptons, a(pp —> Z'X) x BR(Z / —> £ + l~) ) as a function of the three param- 
eters M z >, gy, gx), an d compare it with the limits established by the CDF and DO 
experiments. 

On the theory side, we perform the calculation at NLO in QCD (and LO for the 
EW part), using the NLO MSTW08 PDF sets [27]. In the calculation of the total 
width Y z > we include the following channels: Z 1 — > //, and Zh ) where h is 

the SM Higgs boson and / are the SM fermions of Tab. 1, with the exception of the 
right-handed neutrinos, which we take to be heavier than M^//2. The ratio Tz> /Mz> 
is pretty constant over the whole range of masses of interest, and is of order a few 
percent for GUT-favored Z f couplings, and of course smaller for more weakly coupled 
Z'. 

For the Tevatron experimental limits, we use the most recent available results 
from CDF (on Z' e + e~ [17] and Z' [19]) and DO (on Z' e + e~ [18]). 

They directly provide the 95% CL bounds on the product a(pp — > Z f X) x BR(Z / — ► 
based on 2.5, 2.3, 3.6 fb" 1 of data, with 27 38%, 13 40%, 17 22% total 
acceptances respectively, growing from smaller to larger values of M z >- Notice that, 
although DO data refer to a higher integrated luminosity, the acceptance is smaller, 
making the DO bounds a little weaker than those from CDF. 

When compared with the computed cross-section, the experimental limits pro- 
duce, for each value of M z ^ a 95% CL exclusion region in the (gy,gx) plane, in 
analogy with the EWPT case. 

The other question we want to address is the following. At what combined values 
of center-of-mass (CoM) energy and integrated luminosity may we expect the LHC 
to start having a chance of discovering a Z f (at least of the kind discussed in this 
work), taking into account all the experimental bounds from EWPT and Tevatron 
direct searches? What region of parameter space that has not been already ruled 
out could be accessible for different luminosities and energies in the first LHC runs? 
Considering the fact that Z f signals are among the cleanest and easiest ones in the 
search of new physics, our analysis may also be used as a benchmark point when 
discussing the integrated luminosities that are worth collecting at each energy to 
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actually probe new physics. 

At the moment, the program for the first year of LHC running [4] consists in a 
first run at low energy (y/s = 7TeV) and low luminosity (< 100 pb _1 ), followed by an 
upgrade in energy (y/s < lOTeV), with a collected luminosity up to 200-1-300 pb _1 . 
At such low 7 energies and luminosities, the constraints from Tevatron direct searches 
and EWPT will play a crucial role in identifying the allowed region of parameter 
space that can be probed and the time scale required to have access to it. 

Being the LHC a hadron collider, the region of parameter space accessible to 
the LHC is similar in shape to the corresponding one at the Tevatron. For relatively 
light Z* {Mz> < 800 GeV), since the strongest constraints come from Tevatron direct 
searches, we expect the LHC to turn into a discovery machine as soon as it becomes 
sensitive to regions of parameters not yet excluded by the Tevatron. However, while 
the higher energy is clearly a big advantage for intermediate Z' masses of several 
hundreds GeV, for lighter masses the low luminosity may be a crucial limiting factor 
in the early LHC phase. For heavier Z' masses, such as those relevant for GUT 
models with X = B — L, the discussion is more subtle than in the case of universal 
models discussed in [2]. If the Z' couples significantly to electrons, then generically 
EWPT outperform the Tevatron, and the LHC must wait for higher energies and 
luminosities to become sensitive. If, on the other hand, the Z' couples to muons much 
more than to electrons, then the constraints from EWPT are considerably weaker, 
and the LHC has a significant discovery potential already at a very early phase. 

To turn these considerations into more quantitative statements, we perform a 
basic analysis along the lines of the one described before for extracting the Tevatron 
bounds. In the present case we consider the range y/s — 7 10 TeV for the pp 
CoM energy, luminosities in the range 50 pb _1 1 fb" 1 , and calculate the product 
a(pp -> Z'X) x BR(Z' £+£-) for M z , = 200 3000 GeV, at the same order 
in perturbation theory as in the Tevatron case. At the same level of precision, we 
also compute the SM Drell-Yan (DY) differential cross-section, which constitutes the 
main source of background. 

To gain some approximate understanding of the acceptances for signal and back- 
ground at different values of the invariant mass Mg+i- of the £ + £~ pair, and of the 
possible model-dependence of the former, we performed a simple study, scanning 
over different non-universal minimal models: its conclusions are essentially the same 
as for the universal minimal models of [2]. Since our computed values of the ac- 
ceptance are compatible with those of ref. [28], we adopt their Fig. 2 for the rest 
of our LHC study. More refined studies, however, could take into account also the 
model-dependence of the acceptance, which may not be negligible for Z f searches at 
relatively small masses. 

To estimate the 5a discovery reach of the early phase of the LHC [28, 29, 30], 

7 Of course, with respect to the LHC design parameters. 
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Figure 2: Electrophilic model, X—B—?>L e . Left: The GUT-favored region in the (gY,9x) 
plane and the effect of RGE running from the boundary conditions ajj (10 16 GeV) = 1/24 
and gy(10 16 GeV) = for the two cases (i) and (ii) discussed in the text (higher and lower 
black curves). Right: The shaded regions are allowed at 95% CL by EWPT for M Z ' = 200, 
500, 1000, 1500, 2000, 2500 and 3000 GeV (from inner to outer). 



we compare the events due to a generic Z f signal to the events from the SM-DY 
background in a 3% interval around the relevant values of the dilepton invariant 
mass 8 . We then require the signal events to be at least a 5a fluctuation over the 
expected background, and in any case more than 3 events. This rough statistical 
analysis is enough to get an approximate answer to the questions we want to address. 
We leave a more careful analysis to the experimental collaborations ATLAS and CMS, 
which have control on all the information needed to perform it in an accurate and 
reliable way. 



4. Benchmark models 

4.1 Electrophilic model: X = B — 3L e 

The first model we consider [corresponding to A e = 3, A M = A r = and j3 — 1 in 
Eq. (2.4)] is associated to the gauging of a linear combination of X = B — 3L e and 
the weak hypercharge Y. In this case, when the Z' is 'pure B — 3L e \ it couples only 
to quarks and to first-generation leptons; in particular, the couplings to muons and 
taus only appear when the mixing with Y : parametrized by the effective coupling <?y, 
is non- vanishing. For several aspects this case is similar to the universal case with 
X = B — L, with the exceptions of a stronger coupling to electrons than quarks and 
a suppression of the couplings to muons and taus. The construction of explicit GUT 
models, with a symmetry breaking chain that leads to the present example, would 

8 This conservative assumption is compatible with the expected energy resolution, even in this 
early phase, and with the fact that, for GUT-favored values of the coupling constants, Tz>/Mz> ~ 
few%. 
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Regions allowed by Tevatron Regions allowed by Tevatron 




Figure 3: Electrophilic model, X=B— 3L e . The shaded regions are allowed at 95% CL 
by dielectron (left) and dimuon (right) searches at the Tevatron for M z > = 200, 400, 600, 
800 and 1000 GeV (from inner to outer). 



deserve a separate study, which goes beyond the aim of the present paper. Since 
we are not aware of any fundamental obstruction, we again identify, for illustrative 
purposes, a GUT-favored region in the (gy,gx) plane, performing the RGE evolution 
from Mjj ~ 10 16 down to the weak scale, with boundary conditions varied in the 
interval (2.24). The coefficients of the beta functions in this case are: 
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where: model (i) corresponds to the SM field content plus right-handed neutrinos and 
two scalar Higgs fields, singlets under the SM gauge group, with charges L e — 1 and 
L e = 2; model (ii) corresponds to the MSSM field content plus the scalar superfields 
of three right-handed neutrinos and four Higgs bosons with charges L e = ±1 and 
L e = ±2. The result is plotted on the left-hand side of Fig. 2, which also shows the 
effect of the running (with the generated kinetic mixing) in the case we start from 
the representative values au(Mu) = 1/24 and g Y (Mu) = in both cases (i) and (ii). 
The shape of the preferred region is similar to the X = B — L case. 

The bounds from EWPT can be plotted on the same plane for different choices 
of the Z' mass, see the right-hand side of Fig. 2. Except for a small squeeze in the 
shape of the ellipsoidal curves representing the bounds, induced by the change in the 
relative couplings between electrons and quarks, the bounds are similar to those for 
the universal X = B — L case, since in EWPT most of the observables do not involve 
muons or taus. 

The absence of couplings to muons in Z f x shows up more evidently in the bounds 
from direct searches at the Tevatron, see Fig. 3. Indeed, while the e + e~ channel shows 
a roughly isotropic constraint, independent of the X-Y nature of the Z 7 , as expected, 



-15- 



/7=7TeV, Lum= 100 pb 1 , M z . =0.8 TeV 




gy 

1 /7=7TeV, Lum= 100 pb 1 , M z - = 1.0TeV 




^7=7 TeV, Lum= 100 pb 1 , M z . =0.8 TeV 



/ 

\ 



CDF|K 



\ I 
\ 



-0.4 



-0.2 



0.0 



0.2 



0.4 



By 



0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



/7=7TeV, Lum=100pb-\ M z . = 1.2 TeV 

. . , ^ . . ! . ^ . , . . _ 

X=B~3L e 



-0.4 



-0.2 



feWPT 

o!o 



0.2 



0.4 



Figure 4: Electrophilic model, X—B—3L e . Regions allowed at 95% by precision tests 
(red), by Tevatron e + e~ data (blue), by Tevatron data (green), compared with 

the region not accessible at the LHC (yellow) with = 7 TeV, 100 pb _1 of integrated 
luminosity for Mz> =800, 800, 1000, 1200 GeV respectively. Regions are not shown when 
they would fill the whole plot. The first figure refers to dielectrons searches at the LHC, 
the other three to dimuon searches. The allowed region where the LHC can make a 5a 
discovery is textured. 



the channel shows no constraints for a Z' that is pure B — 3L e . Notice that, on 

the contrary, when B — 3L e mixes with Y ) the channel can produce stronger 

bounds than the e + e~ one; this effect is due to the fact (see Table charges) that the 
B — 3L e and Y contributions to Qz> for the electrons tend to cancel each other in 
the region g Y ~ — (3 6)gx (the lower spikes on the left-hand side of Fig. 3), this 
cancellation being absent for muons that have only the Y contribution to Qz f - 

The shapes of the regions accessible to the LHC are identical to those excluded 
by the Tevatron, only the size is different, depending on the available luminosity 
and CoM energy. The results for the early phase are similar to the B — L case: 
Fig. 4 shows the accessible regions for different masses and channels after 100 pb _1 
of integrated luminosity at 7 TeV, while Fig. 5 refers to 200 pb _1 at 10 TeV. As for 
the B — L case, with the early low-energy low-luminosity run, only very tiny regions 
of parameter space are available for discovery. The situation improves a bit only in 
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Figure 5: Electrophilic model, X=B—3L e . As in the previous figure, but for a better LHC 
with yfs = 10 TeV and 200 pb -1 of integrated luminosity. The left (right) columns refers 
to e + e~ (/i + n~) searches. The Z' mass in the different rows is 600, 900, 1400, 1700 GeV 
respectively. The allowed region where the LHC can make a 5a discovery is textured. 



the region of couplings smaller than those preferred by unification, after the first runs 
at higher energies. An interesting feature is that, at variance with flavor-universal 
models such as the one based on X = B — L, the two channels probe complementary 
regions of parameter space. 
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Figure 6: Electrophilic model, X=B— 3L e . We assume Mz> = 240 GeV, such that the 
CDF e + e~ excess can be reproduced within the blue region at 95% CL. The other shaded 
regions are allowed at the same CL from precision tests (red), CDF dimuon searches (green), 
DO dielectron searches (yellow, hiding behind the blue region). The intersection gives the 
common allowed region: it is shown in textured gray and is non-empty; hence the CDF 
excess is compatible with all bounds. 

The main peculiarity of the electrophilic model with X = B — 3L e is the asymme- 
try between the electron and muon channels at hadronic colliders. This asymmetry 
is able to produce sharp signals in one channel while hiding completely in the other. 
We illustrate this interesting feature using the recent excess [17] observed by the 
CDF collaboration in the dielectron spectrum, at an invariant mass near 240 GeV, 
which however was not observed in the dimuon channel [19]. The question we want 
to answer here is the following: can such an excess be explained with a narrow Z f 
resonance, compatibly with the negative results from the dimuon channel and with 
the constraints from EWPT and from flavor physics? 

The excess would correspond to a resonance with mass around 240 GeV and 
despite the 2.5a significance (which amounts to more than 3a considering only that 
single bin), the coupling to the SM fermions should be small {gz> < 10 -2 ). Since the 
candidate resonance mass is known, we can plot all the different bounds from EWPT 
and Tevatron searches together in the (gy^x) plane: the result is shown in Fig. 6, 
which displays the allowed regions from EWPT (red), CDF e + e~ (blue) and 
(green) searches, DO e + e~ (yellow) searches [18]. In view of the 2.5a excess, the CDF 
e + e~ region does not include the origin (i.e. the SM limit). Notice also that there is a 
non-trivial common region of couplings (textured in gray) which is not ruled out by 
any experiment and is compatible with the observed excess. For such a low-energy 
resonance, the first LHC runs will not be very competitive with respect to the higher 
luminosity data collected at the Tevatron. The values of the allowed couplings are 
quite smaller than those preferred by unification, which suggests that the would- 
be Z f will hardly unify with the other SM gauge coupling at higher energies. It is 
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GUT-favored region Regions allowed by EWPT 




Figure 7: Muonphilic model, X=B—3L /Ji . Left: GUT-favored region in the (gx^gy) plane 
and the effect from the RGE running with boundary conditions ajj (10 16 GeV) = 1/24 and 
gy(10 16 GeV) = for the two cases (i) and (ii) discussed in the text (higher and lower 
black curves). Right: The allowed regions at 95% CL from EWPT for M z > = 200, 500, 
1000, 1500, 2000, 2500 and 3000 GeV (from inner to outer). 

however remarkable that it is possible to fit such an 'exotic' flavor-breaking excess 
with a quite simple, minimal, renormalizable and FCNC-free new physics model. 

4.2 Muonphilic model: X = B — 3L^ 

The second type of Z 1 model we consider [corresponding to A M = 3, A e = A r = and 
(3 — 1 in Eq. (2.4)] is associated to the gauging of a linear combination of X = B—ZL^ 
and Y. This model shares many properties with the previous one for what is related 
to unification and to the signatures at hadron colliders (after the obvious exchange 
of the electron and muon channels). On the other hand, the possibility of having 
a suppression or even the complete cancellation of the couplings to electrons alters 
completely the bounds from EWPT, making them very weak. As we will see in a 
moment, this feature may turn out to be crucial for the early LHC studies. 

For what concerns unification, the story is exactly the same as in the previous 
case: the beta functions are indeed identical, since we just exchanged the role of the 
first and second lepton generation. The plot of the favored region in the (gy,gx) 
plane is thus the same as before and is displayed on the left-hand side of Fig. 7. On 
the right-hand side of the same figure, also the bounds from EWPT for different Z' 
masses are shown. The plot is now completely different from those for X = B — L and 
X = B — 3L e . First of all, there is almost no bound on the Z' masses along the g Y = 
axis: these are indeed pure B — 3L M models, which do not couple to electrons (thus 
no bounds from LEP2) and do not mix with the SM Z boson (thus no constraints 
from LEP1 and APV). The only constraints come from DIS experiments, such as 
NuTeV, and from (g — 2) M , but the bounds are very weak. Notice however that pure 
B — 3L M models are not stable under RGE and in general some mixing with the Y 
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Figure 8: Muonphilic model, X—B—ZL^. Allowed regions at 95% CL from dilepton (left) 
and dimuon (right) searches at Tevatron for M Z / = 200, 400, 600, 800 and 1000 GeV (from 
inner to outer). 



will be produced; still, the region around gy = is weakly constrained by EWPT, 
which makes hadronic colliders more sensitive to this particular type of models. 

The bounds from the Tevatron experiments are reported in Fig. 8 and are pretty 
similar to those of the previous section, with dielectron and dimuon channels ex- 
changed. Notice that now the muon channel is the most powerful for probing the 
]]y — region, where EWPT are weaker. Again, the plots for the LHC share the 
same shapes as those from the Tevatron, the size being different and dependent on 
the energies and luminosities considered. 

The attractive feature of this model is thus evident. We have seen that there 
is a region of couplings where EWPT are particularly weak (g Y — 0) and where 
we may expect the dimuon searches at hadron colliders to be very powerful instead. 
We also know that for resonances above 1 TeV the higher energy of the LHC can 
easily outperform the Tevatron searches and bounds, even after collecting a very small 
integrated luminosity. There are then regions allowed by EWPT and Tevatron where 
the LHC can make a discovery after collecting only few pb _1 of data. This possibility 
is illustrated in Fig. 9, which shows the combination of bounds from EWPT, Tevatron 
bounds and the LHC reach for M z > = 800, 1000 and 1200 GeV. The textured regions 
correspond to Z' models not excluded by existing experimental bounds and accessible 
for a 5a discovery, already in the first phase of LHC: as representative values, we 
take 50 pb _1 of integrated luminosity at yfs = 7 TeV. Notice that for Mz> = 1 TeV 
the accessible region would actually include the one preferred by unification (even 
considering the mixing induced by the RGE running) . In the language of [20] , these 
models can be considered to be supermodels. This very simple possibility was not 
considered in [20] because of possible issues with FCNC, but in Sect. 2.2 we showed 
that such flavor problem is automatically solved in our class of models. 
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Figure 9: Muonphilic model, X=B—3L fjL . Regions allowed at 95% C.L. by precision 
tests (red), Tevatron e + e~ data, Tevatron data (green) compared to the region not 

accessible at 5a level by the LHC dimuon searches (yellow) with y/s = 7 TeV, 50 pb _1 of 
integrated luminosity. We assume M z > =800 (upper plot), 1000 (medium) and 1200 GeV 
(lower). The allowed region where the LHC can discover a 5a signal is textured. 



4.3 Hadrophobic model: X = — L T 

Finally, the B — 3L T model does not lead to clean peaks in the lepton invariant masses 
at hadron colliders. It is weakly constrained by hadronic Tevatron data and by the 
non-observation of matter effects in <-> v T atmospheric oscillations. Therefore we 
prefer to consider a different last example [corresponding to A r = — A M = 1, A e = 
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Figure 10: Hadrophobic model, X=L M — L T . Left: GUT-favored region in the (gx,9Y) 
plane and the RGE running from a^(10 16 GeV) = 1/24 and gy (10 16 GeV) = for the two 
models (i) and (ii) discussed in the text (higher and lower black curves). Right: Regions 
allowed at 95% CL by precision tests for M z > = 200, 500, 1000, 1500, 2000, 2500 and 
3000 GeV (from inner to outer). 

and — in Eq. (2.4)], where the Z f couples to a linear combination of Y and 
X = L M — L T . The peculiarity of this model is that, in the pure L M — L T phase, the 
Z' couples neither to electrons nor to quarks, hiding completely from EWPT and 
direct collider searches. In this particular case the Z' can be extremely light, the 
only bounds coming from the muon magnetic moment. When the mixing with Y is 
considered, the Z' couplings to quarks and electrons switch on and the usual bounds 
apply again, albeit suppressed by the mixing angle 9' when the latter is small. Unlike 
the previous cases, we thus expect a dark region along the g Y = axis, blind both 
to bounds and to discoveries at colliders. 

Also the properties under RGE running are different in this model. Indeed, no 
kinetic mixing is generated at one loop if it vanishes at tree level: 
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The absence of kinetic mixing makes the axis gy — an attractor under the RGE 
evolution, which means that pure L M — L T models are stable under quantum effects 
and do not mix appreciably with F, even in the case when small threshold corrections 
are present at the GUT scale. This can be seen by looking at the RGE for g Y near 
9y = 0: 

It = i& ( 2 w * + bxxgl + ° {q2y) ) * 

which shows that dgyjdt vanishes at gy = and has the same sign as gy, so that 
the latter flows to zero at low energies. The GUT-favored region in this case is 
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Figure 11: Hadrophobic model, X=L M — L T . Regions allowed at 95% CL by dilepton 
(left) and dimuon (right) searches at Tevatron for M z > = 200, 400, 600, 800 and 1000 GeV 
(from inner to outer). 
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Figure 12: Hadrophobic model, X—L^—L T . Comparison of 95% CL bounds from EWPT 
(red) and Tevatron (blue for e + e~ and green for and discovery reach at the LHC 

(yellow) with = 10 TeV, 200 pb _1 of integrated luminosity for Mz> =1 TeV for dielec- 
tron and dimuon searches respectively. The textured region is the one accessible to a 5a 
discovery and compatible with existing bounds. 

represented in Fig. 10, the absence of kinetic mixing makes the region symmetric 
around the vertical axis. 

The constraints from EWPT at various masses are displayed in the same Fig. 10. 
As expected, no useful constraints are present for a Z' associated to pure L M — L T . 
Away from the vertical axis the constraints are basically on the 'y-amount' gy of 
the Z', almost independently of gx- 

The bounds from the Tevatron searches are plotted in Fig. 11. In the dielectron 
channel, the region around the unconstrained g Y = line (pure L M — L T ) becomes 
larger with increasing g X) because larger g x increases the couplings to muons and 
taus, thus dumping the branching ratio to electrons. 
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The dimuon bounds are not symmetric under g Y —> —gy, and are less con- 
straining for gx ~ 9y, because of a partial cancellation between the Y and L M — L T 
components of the Z' couplings to muons. 

Coming to the LHC discovery reach, there are regions accessible in the early 
phase, both in the e + e~ and channels. These are illustrated in Fig. 12, where 

we plot, as a representative example, the accessible regions for a 5a discovery at the 
LHC after 200 pb" 1 at y/s = 10 TeV, for M z < = 1 TeV. The muon channel seems 
quite powerful, exploring a sizable region of the (gy,gx) plane for these masses, 
except of course the dark region around gy — 0. For smaller masses the Tevatron 
bounds start dominating, while at larger masses the bounds from EWPT eventually 
overcome the early LHC reach. 

5. Conclusions 

Starting from the SM augmented by three right-handed neutrinos, to account for neu- 
trino masses, we explored the most general extra U(l) that can be gauged compatibly 
with the following requirements: the theory must be renormalizable, anomaly-free, 
compatible with the observed fermion masses and mixings and with the absence of 
extra violations of flavor, beyond the SM CKM matrix and those induced in the 
lepton sector by the light neutrino masses. The most general Z' boson of this kind 
was found to be a linear combination of the hypercharge F, B — 3L e , B — 3L M , 
B — 3L r . Although this statement is not new, we showed that it still holds now 
that neutrino masses and mixings have been observed. Unlike in the quark sector, 
where individual baryon flavors cannot be gauged, this is still possible in the lepton 
sector, as the resulting flavor-changing neutral currents are suppressed by the small 
neutrino masses, as in the SM. 

A new Z' gives clean easy signals at the LHC but also in previous experiments, 
thereby we computed the allowed parameter space comparing the present bounds 
from EWPT and the Tevatron to the LHC capabilities. It is important to keep in 
mind that a given Z 7 , associated with a definite direction in the four-dimensional 
space defined by Y and B — 3L a (a = e,/x,r), is specified by two parameters: its 
mass M z > and its effective gauge coupling g z >. The often-reported limits on the Z' 
mass [1] tacitly assume some specific arbitrary value of g z > . Although we identified 
the ranges of g z > suggested by grand unification, making use of the appropriate RGE 
from the GUT scale to the weak scale, this is just a theoretical hypothesis. Both 
parameters must be kept in a complete phenomenological analysis, and both are 
important. The left-hand side of Fig. 13 shows a typical case: the region where the 
early LHC can outperform all previous experiments corresponds to M z > > 500 GeV 
and g Z t significantly smaller than the GUT-favored values. Indeed, if the Z' is light, 
y/s is not important, and the early LHC cannot compete with the ~ 50 times larger 
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Figure 13: Early LHC 5a discovery reach compared with bounds from EWPT and Teva- 
tron. Left: pure B — 3L e . Right: pure B — 3L M . For values of the LHC energy and 
luminosity higher than those indicated in the figure, the sensitivity in gx roughly scales as 
(luminosity) -1 / 2 , while the sensitivity in Mz' roughly scales as y/s. 

luminosity already accumulated at the Tevatron. Furthermore, too large values of 
gz' /M z > have already been excluded indirectly by precision data. 

We focussed on the following representative possibilities. 

1. The elect rophilic B — 3L e , possibly mixed with the hyper charge Y. The left- 
hand side of Fig. 13 compares the LHC sensitivity with the bounds from Teva- 
tron and from EWPT, it has the typical shape of similar graphs for generic 
Z 1 models. The main novelty with respect to the B — L case studied in [2] 
is the possibility to fit the CDF excess in the e + e~ channel with M z > ~ 240 
GeV compatibly with all other constraints, in the region of parameter space 
shown in Fig. 6. Being light, such Z 1 will be tested at the Tevatron with more 
luminosity sooner than at the LHC. 

2. The muonphilic B — 3L M , possibly mixed with the hypercharge Y. The main 
novelty is that this possibility is very weakly constrained by precision data when 
the mixing with Y is small: the two experiments most sensitive to the pure 
B — 3L M case, g — 2 and NuTeV, show anomalies, but it is not easy to fit them 
and the impact on the allowed parameter space is negligible. As illustrated on 
the right-hand side of Fig. 13, in the case of pure B — 3L M , a significant region 
of parameter space can therefore be explored even by the initial stage of the 
LHC with reduced energy and luminosity. 

3. Finally, the hadrophobic L M — L r , possibly mixed with the hypercharge Y, can 
give detectable signals at the LHC, but not easily in the early phase. Unlike 
the first two cases, mixing with hypercharge is not automatically generated by 
the RGE running. 
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